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We propose a method to extract the projected power spectrum of density perturbations from the 
distortions in the cosmic microwave background (CMB). The distortions are imprinted onto the 
CMB by the gravitational lensing effect and can be extracted using a combination of products of 
CMB derivatives. We show that future CMB experiments such as Planck will be able to extract 
the power spectrum with high statistical significance over two orders of magnitude in angle. The 
method proposed here traces dark matter directly to higher redshift (up toz~ 1100) and larger scale 
(few Gpc) than any other currently known method. It also traces large scale structure in the linear 
regime, allowing simple interpretation in terms of cosmological models. By providing additional and 
complementary information to the one from the primary CMB analysis it will strengthen further 
the scientific return of future CMB experiments. 
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The standard paradigm of modern cosmology is the ex- 
panding universe model, where the universe started with 
a hot big bang and has been expanding ever since. The 
universe is predicted to be homogeneous and isotropic 
on large scales with initially small fluctuations superim- 
posed on that background. These fluctuation have been 
amplified by gravity to become the present day large scale 
structure. This paradigm, while successful in explaining 
all the observations, still leaves many questions unan- 
swered. Among these are the value of cosmological pa- 
rameters, such as the density of the various components, 
the nature of dark matter and the power spectrum of 
initial perturbations. 

Observational efforts are necessary to provide the an- 
swers to the above questions. It has long been recog- 
nized that the power spectrum of dark matter fluctu- 
ations would provide many of these answers, as it de- 
pends both on the initial power spectrum as well as on 
the rate of growth of perturbations, which are determined 
by the cosmological model and the nature of dark matter 
(e.g. hot or cold). Its measurement has been the pri- 
mary target of many observational surveys in the past 
and this trend will also continue in the future. Most 
of these have used galaxy positions and redshifts to map 
the universe. The ongoing surveys such as SDSS and 2dF 
hope to measure the power spectrum accurately on scales 
up to several hundred megaparsecs. However, there are a 
number of possible uncertainties that may make this goal 
more difficult than expected. In particular, the fact that 
galaxy formation is such a poorly understood process 
makes the connection between mass and light uncertain. 
The relation between dark matter and galaxy density 
perturbation could be nonlinear, scale-dependent, non- 
local or stochastic even on very large scales . Another 
complication is the fact that except on very large scales, 
the surveys are measuring the clustering in the nonlinear 



regime, which makes the interpretation and statistical 
analysis of the data more complicated. 

For these reasons other, more direct, tracers of dark 
matter have been investigated, most notably velocity 
flows ||, clustering of Ly-a forest [|| and weak lensing 
jlj . Each of these has uncertainties associated with them. 
Weak lensing, although not yet detected, offers perhaps 
the best chance of measuring accurately the power spec- 
trum of dark matter over a large range of scales. By 
measuring the distortions in the shapes of background 
galaxies deformed by the mass distribution along the line 
of sight one can determine the clustering fluctuations on 
scales from 1 up to 100 7i _1 Mpc Weak lensing has 

the advantage of tracing the dark matter directly and 
therefore avoids the uncertainties connected with galaxy 
clustering. Still, there are uncertainties associated with 
weak lensing as well, the most important being the poorly 
known redshift distribution of background galaxies and 
their possible intrinsic alignment. 

In this letter we propose a new method to measure the 
projected dark matter power spectrum by using the weak 
lensing distortions of the CMB |t]]. We show that grav- 
itational lensing induces characteristic distortions in the 
pattern of the CMB field, which allows one to reconstruct 
the projected dark matter density field Q. The recon- 
structed projected mass density can be used to measure 
the mass power spectrum by averaging over independent 
CMB patches. The method proposed here offers several 
advantages over the other methods discussed above. It 
traces the dark matter directly and does not depend on 
the assumptions of how light traces mass. Unlike weak 
lensing of galaxies it does not suffer from the possible 
intrinsic alignment of background galaxies or their un- 
certain redshift distribution. It is able to recover the 
power spectrum over a large range of scales, roughly be- 
tween 10 Mpc and 1 Gpc. As such it is sensitive to scales 
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larger than any other survey, including the SDSS and 
2dF. Furthermore, the observed field is a projection of 
density up to 2 ~ 1100 and one is therefore tracing dark 
matter to higher redshifts than with any other method. 
The last two features also imply that with the proposed 
method one is tracing the perturbations predominantly in 
the linear regime and no nonlinear corrections arc neces- 
sary down to the smallest scales that are still observable. 
In many ways the method proposed here offers the same 
advantages that make the CMB anisotropics such a pow- 
erful test of cosmological models, except that it is pro- 
viding complementary information on the projected dark 
matter power spectrum. Clearly, combining the standard 
CMB analysis with the present one will further enhance 
the power of CMB to constrain the cosmological mod- 
els || and provide additional motivation for future CMB 
experiments. 

To provide a clear physical interpretation we will only 
discuss the small scale limit of the method in this letter 
(see H|l0| for all-sky generalization). The observed CMB 
temperature in the direction 9 is T{9) and equals the (un- 
observable) temperature at the last scattering surface in 
a different direction, T(9 + 59), where 59 is the angu- 
lar excursion of the photon as it propagates from the 
last scattering surface to us. In terms of Fourier com- 
ponents we have T{9) = (2n)- 2 J d 2 l e u < e+50 ^ f{l).To 
extract the information on the deflection field 50 we con- 
sider derivatives of the CMB temperature. If the CMB 
is an isotropic and homogeneous Gaussian random field 
then different partial derivatives are statistically equiva- 
lent and their spatial properties are independent of po- 
sition. Lensing will distort these two properties of the 
derivatives. The derivatives of the temperature field are 
to lowest order, 



dT 



59) = (5 ab + $ ab )f b (0 + 50), (1) 



where $> ab = is the symmetric shear tensor. The 
components of the shear tensor can be written in terms 
of the convergence n = — (& a a + ®bb)/2 and the two shear 
components 71 = -($ oa - $b&)/2 and 72 = -$ a &. 

Next we consider the quadratic combinations of the 
derivatives and express them in terms of the unlensed 
field to lowest order in the shear tensor, 
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from above (S) = 2k, (Q) = 2 7 i, (U) = 2 72 . Physically, 
k stretches the image and makes the derivatives smaller. 
Similarly, shear produces anisotropy in the derivatives 
and can be extracted by considering the particular com- 
bination of derivatives defined above. From Q and IA we 
can form two rotationally invariant quantities in Fourier 
space 
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where as was defined so that in the absence of lensing 
(<S) = (Q) = (U) = 0. In the presence of lensing it follows 



£{l) = e(I)cos(20i) 
B{1) = Q(i) sin(2&) - 



-U(l)sm(24>i) 
■U(l)cos(24>i). 



(3) 



The average of the scalar field £ is {£) = 2k. The average 
of the pseudo-scalar field (£>) vanishes in the large scale 
limit because gravitational potential from which shear is 
generated is invariant under the parity transformation. 
The convergence k can be reconstructed in two indepen- 
dent ways, either from S oi £, while B can serve as a 
check for possible systematics. 

The convergence k can be written as a projection 
of gravitational potential k = / Q Xo g(x, Xo)V 2 0(x)dx 
f|,||. Here xo is the comoving radial coordinate at 
recombination and g is the radial window, defined as 
9(XiXo) — T ^r(xo)~ X ' > ' H ere r (x) is the comoving an- 
gular diameter distance, defined as if -1 / 2 sin K x l 2 x, X> 
(- J ftT)- 1 / 2 sinh(- J fs:) 1 / 2 x for K > 0, K = 0, K < 0, re- 
spectively. The curvature K can be expressed using the 
present density parameter f2o and the present Hubble 
parameter Hq as K = (flo — 1)Hq. In general Slo con- 
sists both of matter contribution Q m and cosmological 
constant term 

The angular power spectrum of convergence Cf K has 
ensemble average H 



X " 9 2 (x,Xo) 
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where a(x) is the expansion factor and we used Poisson's 
equation to express the convergence directly in terms of 
time dependent matter power spectrum Pg{k, t). Fig- 
ure ^ shows the contribution of the 3-d power spectrum 
Ps(k) to Cf K as a function of k for different I. It is 
a relatively broad function of k. Different cosmological 
models give somewhat different correspondence between 
I and k, but in general I — 10 probes scales around 
A = 2n/k = l/i _1 Gpc, while I = 1000 contribution 
probes scales around A = 30/i _1 Mpc. We will see below 
that this spans the observable range of Planck satellite. 

From S and £ one can form three different angular 
power spectra, Cf s , Cf £ and Cf £ , defined as C ; ww ' = 
[W(i)*W(i') + W(l')*W(l)]/2 5 lv . Their expectation 
value is (C z ww ') = 4^Cf K + N™ w . Here N™ w are 
the noise power spectra for S, £ or their cross-term, while 
Wi is the window function describing the signal degra- 
dation because of finite angular resolution and detector 
noise ls|. The window changes smoothly from unity at 
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FIG. 1. Logarithmic contribution to Cf " as a function of 
k for I — 10, 100, 1000 (the normalization is arbitrary). The 
models are flat CDM model (dotted), open CDM model with 
f2m = 0.3 (dashed) and cosmological constant model with 
with fl m — 0.3 (solid). All the models have F = fl m h = 0.21. 



low I to zero at high I. For MAP the transition occurs at 
I ~ 300, while for Planck the window is close to unity up 
to I ~ 1000. The noise power spectra have contribution 
from intrinsic CMB fluctuations and detector noise. In 
the absence of detector noise it is the correlation length 
of the CMB £ that governs the amplitude of noise on 
large scales. Each patch of the sky of size £ 2 contributes 
one independent measurement with a variance of order 
unity. By averaging over many independent patches we 
can reduce the noise on large scales. Detector noise and 
beam smoothing degrade this and the CMB field has to 
be smoothed at the angular scale where the detector noise 
exceeds the CMB signal S . This increases the CMB cor- 
relation length and the amplitude of noise power spec- 
trum. Beam and noise characteristics of Planck (0.12° 
FWHM and N TT = (O.Oljulf) 2 ) make the influence of de- 
tector noise negligible, as there is very little CMB power 
on small scales which are not accessible by Planck. The 
noise spectra for SS, ££, BB and S£ are shown in top of 
figure |. On large scales N ££ w N m w N ss /2 > N S£ . 
Also shown is the typical convergence power spectrum 
4Cf K , which is significantly below Planck noise for S 
or £ . This shows that the reconstruction can only be 
achieved in a statistical sense by averaging over indepen- 
dent multipole moments. For MAP with 0.21° FWHM 
and N TT = (0.11/iK) 2 the noise power spectra are a fac- 
tor of 5 higher than for Planck @] . From the measured 
Cj vw we subtract the noise power spectra N^ w and 
are left with an estimate of C; = 4M / /Cf K . 

To combine the three estimates of C\ obtained from 
SS, ££ and S£ correlations in an optimal way, we con- 



sider their covariance matrix. If the CMB noise can be 
considered Gaussian then the covariance matrix for the 
estimated power spectra can be expressed in terms of 
the noise power spectra. We have verified using Monte 
Carlo simulations that the Gaussian approximation is 
an excellent one in the large scale limit H, despite the 
fact that the underlying fields are not Gaussian (they 
are constructed using squares of the temperature field). 
This is a consequence of the central limit theorem, as 
there are many uncorrelated patches that contribute to 
the low I modes. The diagonal terms of the covari- 
ance matrix are given by Cov[(C ww ) 2 ] = ^{N^f 
and Cov[(C S£ ) 2 } = ^[N^Nf* + (Nf £ ) 2 ]. The off- 
diagonal elements can be ignored compared to the di- 
agonal terms at low I because Nf £ <C Nf s ,Nf £ . The 
covariance matrix is then diagonal and we can use simple 
inverse noise variance weighting to find the best combi- 
nation 
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Cov 1 [(C S£ ) 2 } is the variance of Q. 

The reconstructed average multiplied with W^ 1 is 
plotted in figure for the Planck experiment using one 
Monte Carlo realization of the sky. We find that the in- 
put power spectrum is recovered up to I ~ 1000. The 
ratio between the two, shown in bottom of figure ^, is 
consistent with unity over this range, while the corre- 
sponding ratio for B is consistent with 0. There is no 
significant bias in S, £ or B. 

Note that even for Planck noise is larger than the sig- 
nal by a large factor, (S/N)' 1 = N ££ /4WiCf K ~ 10 
around I ~ 20 for this model. Therefore we need 
about 2(S/N)i 2 ~ 200 modes to reach S/N = 1 on 
the power spectrum and this in only possible if I > 
(2/f sky ) 1 / 2 (S/N) l ~ Uf^ 2 , where / sky is the sky frac- 
tion that is being observed. This means that we cannot 
successfully recover the power spectrum for I < 14 with 
S/N > 1, so the small scale analysis in this paper is ad- 
equate. However, the information from low I modes can 
still be useful. There may be models which predict a 
large increase in power on very large (Gpc) scales. For 
such models noise may be below the signal and should 
give a detectable signal. Also, low / modes can be cross- 
correlated with other maps to enhance the signal to noise. 
One example is the CMB itself, where a positive detection 
would be a signature of a time dependent gravitational 
potential {LG|-|^. 

For MAP the reconstruction above fails to give a 
positive detection. A more detailed statistical analysis 
is necessary to asses the signal to noise when all the 
information is combined. To asses the overall signal 
to noise we combine the information from all the esti- 
mates X = such that S/N = (X)/(X 2 ) 1 / 2 = 
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FIG. 2. Top: power spectra for noise N? s (short dashed), 
Nf £ (dotted), Nf s (long-dashed) and TVf 8 (dash-dotted) for 
Planck, using cosmological constant model with Q. m = 0.3 and 
(jg = 1. Also shown is power spectrum of convergence 4Cf K 
(solid) for the same model, together with its reconstruction 
from a Monte Carlo simulation. Bottom: ratio of output to 
input power spectrum is plotted for the averaged C; reK recon- 
struction and for the Cf B reconstruction. Also shown is the 
window Wi (dashed) for Planck satellite characteristics. 



approximation may still be valid on large scales because 
of the central limit theorem. The accuracy of subtraction 
is more critical for C ss and C ££ than for C S£ , where the 
noise is lower than the signal for I < 300 (figure ||) . If only 
cross-correlation is used to estimate C KK then the signal 
to noise is becomes roughly 2 /3 of our previous estimates 
and would still be clearly detectable with Planck. Finally, 
there are many consistency checks that can be applied to 
the reconstruction, such as comparing the three estima- 
tors which should all be consistent with each other and 
verifying that C is consistent with 0. For low I one 
may only use N BB to subtract the noise directly, as in 
this limit we have N BB = N ££ = N ss /2. To summarize, 
future CMB experiments have the potential of providing 
another important map of the universe, allowing one to 
trace the large scale distribution of dark matter to higher 
redshifts and larger scales than any other method. 
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by ai = Gij u 2 c [hOl . The resulting signal to noise is 



Y^i otiCi/i^n af On, Y' 2 is maximized. This is achieved 
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For MAP this gives S/N = 3 for fl m = 0.3, cr 8 = 1 model 
with cosmological constant. Other viable models give 
comparable numbers. MAP detection will therefore only 
be marginal, unless the power spectrum on large scales 
turns out to be much larger than expected. For Planck 
the viable models give S/N between 15-25, confirming 
the conclusion above that Planck will be able to extract 
the power spectrum with high statistical significance. 

We have also investigated possible sources of system- 
atics. One is the fact that we have to subtract the CMB 
noise, which depends on the detailed knowledge of the 
unlensed CMB power spectrum, while only the lensed 
power spectrum is known with finite accuracy. We have 
found both effects to be below 0.5% of the noise am- 
plitude of Planck, so they will not significantly affect the 
power spectrum reconstruction for Planck, but will be an 
additional source of error for MAP || . Another possible 
uncertainty is the assumption of the Gaussian statistics, 
which allows us to estimate the noise and subtract it from 
the measurements. It should be noted that even if there is 
a non Gaussian component in the CMB (for example due 
to foregrounds and secondary processes), the Gaussian 
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